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CHARM MIXING IN THE STANDARD MODEL AND BEYOND 



The motivation most often cited in searches for — Z)'^ mixing and CP-violation in 
charm system lies with the possibility of observing a signal from New Physics which 
dominates that from the Standard Model. We review recent theoretical predictions and 
experimental constraints on D^ — D'^ mixing parameters, concentrating on possible effects 
of New Physics. 

PACS numbers: 11.25.Hf, 123.1K 

1. Introduction 

Quantum mechanical meson-antimeson oscillations are sensitive to heavy degrees of 
freedom which propagate in the underlying mixing amplitudes. Comparing observed 
meson mixing with predictions of the Standard Model (SM), modern experimen- 
tal studies would be able to constrain models of New Physics (NP). Yet, extensive 
precission data from B-factories and the Tevatron collider show that the large SM 
mixing succesfuUy describes all available experimental data in and oscilla- 
tions. The only flavor oscillation not yet observed is that of the charmed meson 

, where SM mixing is very small and the NP component can stand out. This 
situation is an exact opposite to what happens in the B system, where S*^ — B^ 
mixing measurements are used to constrain top quark couplings. 

Together with the one loop Standard Model effectJU NP effects can contribute 
to AC = 1 (decays) or AC = 2 (mixing) transitions. In the case of D° — mixing 
these operators generate contributions to the effective operators that change _D" 
state into Z?" state, leading to the mass eigenstates 



where the complex parameters p and q are obtained from diagonalizing the Z?*^ — 
mass matrix. Note that |pp -|- |gp = 1. If CP-violation in mixing is neglected, p 
becomes equal to q, so \Di^2) become CP eigenstates, CP\D±) = ±\D±). The mass 
and width splittings between these eigenstates are given by 
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\Di^2)=p\D')±q\D'^) 



(1) 
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It is known experimentally that D'^ — mixing proceeds extremely slowly, which 
in the Standard Model is usually attributed to the absence of superheavy quarks 
destroying GIM cancellations. 



2. Phenomenology of D° — D° mixing 

Semileptonic decays. The most natural way to search for charm mixing is to 
employ semileptonic decays. It is also not the most sensivite way, as it is only 
sensitive to Rd — {x^ + ?/^)/2, a quadratic function of x and y. Use of the D° 
semileptonic decays for the mixing search involves the measurement of the time- 
dependent or time-integrated rate for the wrong-sign (WS) decays of D, where 
c — > (c via mixing) si~v, relative to the right-sign (RS) decay rate, c —^ st^v. 

Nonleptonic decays to non-CP eigenstates. Currently the most strin- 
gent limits on the mixing parameters arise from the decay time dependent 

—)■ K^TT^ measurements. Time-dependent studies allow separation of the di- 
rect doubly-Cabbibo suppressed (DCS) K^tt~ amplitude from the mixing 
contribution Z?" K^tt^ , 

Tdo^k+tt- e-'^^'lAx-^+f R + \^R,n{y' cos (t)^ x' sin (t))rt + R^Rjjinf (3) 

where R is the ratio of DCS and Cabibbo favored (CF) decay rates. Since x and y 
are small, the best constraint comes from the linear terms in t that are also linear in 
X and y. A direct extraction of x and y from Eq. ([3]) is not possible due to unknown 
relative strong phase Sd of DCS and CF amplitudes, as x' = xcosSD+ysinSD, y' — 
y cos Sd — xsihSd- This phase can be measured independently. The corresponding 
formula can also be writterPlfor D° decay with x' ~x' and Rm ^ ■ 

Nonleptonic decays to CP eigenstates. 13° mixing can be measured by 
comparing the lifetimes extracted from the analysis of D decays into the CP-even 
and CP-odd final states. In practice, lifetime measured in D decays into CP-even 
final state fcp, such as K~ ,tt'^tt~ , cj)Ks, etc., is compared to the one obtained 
from a measurement of decays to a non-CP eigenstate, such as K^t:^ . This analysis 
is also sensitive to a linear function of y via 



t(D ^ K-7r+) 



Rm — 1 



(4) 



where is a CP-violating phase. In the limit of vanishing CP violation ycp — V- 
This measurement requires precise determination of lifetimes. 

Quantum-correlated final states. The construction of tau-charm factories 
introduces new time-independent methods that are sensitive to a linear function of 
y. One can use the fact that heavy meson pairs produced in the decays of heavy 
quarkonium resonances have the useful property that the two mesons are in the 
CP-correlated states^. For instance, by tagging one of the mesons as a CP eigen- 
state, a lifetime difference may be determined by measuring the leptonic branching 
ratio of the other meson. The final states reachable by neutral charmed mesons 
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Standard Model mixing predictions 




Reference Index 



Fig. 1. Standard Model predictions for \x\ (open triangles) and \y\ (open squares). 



are determined by a set of selection rules according to the initial virtual photon 
quantum numbers J^^ — 1 Since we know whether this D{k2) state is tagged 
as a (CP-eigenstate) D± from the decay of D{ki) to a final state Sa of definite 
CP-parity cr = ±, we can easily determine y in terms of the semileptonic branching 
ratios of D±, which we denote Bj^. Neglecting small CP-violating effects, 



^\Bl{V) BXiV), 

A more sophisticated version of this formula as well as studies of feasibility of this 
method can be found in Ref. [3]. 

The current experimental bounds on y and x are SI 

y < 0.008 ± 0.005 , x < 0.029 (95% C.L.) . 



3. Charm mixing predictions in the Standard Model 

The current experimental upper bounds on x and y are on the order of a few times 
10"'^, and are expected to improve in the coming years. To regard a future discovery 
of nonzero cc or ?/ as a signal for new physics, we would need high confidence that 
the Standard Model predictions lie well below the present limits. As was shown in 
[5], in the Standard Model, x and y are generated only at second order in SU(3)f 
breaking, 

X , y - sin^ Oc x [SU{3) breaking]^ , (6) 

where 6c is the Cabibbo angle. Therefore, predicting the Standard Model values of 
x and y depends crucially on estimating the size of SU(3) p breaking. 

Theoretical predictions of x and y within the Standard Model span several orders 
of magnitude^ (see Fig.[T]). Roughly, there are two approaches, neither of which give 



''Horizontal line references arc tabulated in Table 5 of Ref. [6]. 
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very reliable results because rric is in some sense intermediate between heavy and 
light. The "inclusive" approach is based on the operator product expansion (OPE). 
In the rric ^ A limit, where A is a scale characteristic of the strong interactions, 
AM and AF can be expanded in terms of matrix elements of local operators^. 
Such calculations typically yield x,y < 10~^. The use of the OPE rehes on local 
quark- hadron duality, and on A/mc being small enough to allow a truncation of the 
series after the first few terms. The charm mass may not be large enough for these 
to be good approximations, especially for nonleptonic D decays. An observation of 
y of order 10~^ could be ascribed to a breakdown of the OPE or of duality, but 
such a large value of y is certainly not a generic prediction of OPE analyses. The 
"exclusive" approach sums over intermediate hadronic states, which may be modeled 
or fit to experimental datcP. Since there are cancellations between states within a 
given SU{3) multiplet, one needs to know the contribution of each state with high 
precision. However, the D is not light enough that its decays are dominated by a 
few final states. In the absence of sufficiently precise data on many decay rates and 
on strong phases, one is forced to use some assumptions. While most studies find 
x,y < 10~^, Refs. [8] obtain x and y at the 10~^ level by arguing that SU(3)f 
violation is of order unity. It was also shown that phase space effects alone provide 
enough SU(3)_f violation to induce x,y ^ 10~^ [5]. Large effects in y appear for 
decays close to D threshold, where an analytic expansion in SU(3) f violation is no 
longer possible; a dispersion relation can then be used to show that x would receive 
contributions of similar order of magnitude. 

The above discussion shows that, contrary to B and K systems, theoretical 
calculations of x and y are quite uncertaiJ^, and the values near the current exper- 
imental bounds cannot be ruled out. 



4. New Physics contribution to D° — D° mixing 

As one can see from the previous discussion, mixing in the charm system is very 
small. As it turns out, theoretical predictions of x and y in the Standard Model 
are very uncertain, from a percent to orders of magnitude smaller Thus, New 
Phyiscs (NP) contributions can easily stand out. 

In order to see how NP might affect the mixing amplitude, it is instructive to 
consider off-diagonal terms in the neutral D mass matrix, 

V 2 „ MB-En + ie 

where H^'^'^^^ is the effective |AC| — 1 hamiltonian. 

New Phyiscs in |AC| — 2 interactions. Since all new physics particles are 
much heavier than the Standard Model ones, the most natural place for NP to affect 
mixing amplitudes is in the |AC| — 2 piece, which corresponds to a local interaction 
at the charm quark mass scale. Integrating out NP degrees of freedom at some scale 
A, we are left with an effective Hamiltonian written in the form of series of operators 
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New Physics mixing predictions 
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Fig. 2. New Physics predictions for 



of increasing dimension^. Realizing this, it is not hard to write the complete basis 
of those effective operators, which most conveniently can be done in terms of left- 
and right-handed quark fields, 

1=1 

where Ci are the Wilson coefficients, and Qi are the effective operators, 

Q2 ^UH'Jfj^CRULj^CL, Qe =UR"ffj.CRURj^CR, (9) 

Q3 = UlCrUrCl, Q7 = UlCrUlCr, 

Q4 ^ URCLURCL, Qs = UL(Tfj.„CRUL<J^'' CR, 

These operators exhaust the list of possible contributions to AC = 2 transitions. 
Since these operators are generated at the scale /i = A (at which new physics is 
integrated out), a non-trivial operator mixing can occur if we take into account 
renormalization group running of these operators between /i = A and ~ rtic 
scales. This running can be accounted for by solving RG equations obeyed by the 
Wilson coefficient functions, 

^—C{^^)^^^{^^)C{^.), (lo) 

CtlOg/Lt 

where ^"^{n) represents the matrix of anomalous dimensions of operators of Eq. 
A prediction for a mixing parameter x in a particular model of new physics is then 
obtained by computing d (A) for a set of Qi (A) generated by a given model, running 
the RG equations of Eqjf fTU)) and computing matrix elements {D'^\Qi{mc)\D^) . As 
can be seen from Fig. predictions for x vary by orders of magnitude for different 



''Horizontal line references are tabulated in Table 5 of Ref. [B]. 
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models. It is interesting to note that some models require large signals in the charm 
system if mixing and FCNCs in the strange and beauty systems are to be small 
(e.g. the SUSY alignment model). 

New Phyiscs in |AC| — 1 interactions. The local |AC| = 2 interaction 
cannot, however, affect AFd because it does not have an absorptive part. Thus, 
naively, NP cannot affect lifetime difference y. This is, however, not quite cor- 
rect. Consider a D'^ decay amplitude which includes a small NP contribution, 
A[D^ — > n] = A^^'' +A^^\ Here, ^i^^"* is assumed to be smaller than the current 
experimental uncertainties on those decay rates. Then it is a good approximation 
to write y as 

y - E f^Mr^^^n""^ + 2 E Y^^r^^^n""^ ■ (11) 
n n 

The SM contribution to y is known to vanish in the limit of exact flavor SU{3). 
Moreover, the first order correction is also absent, so the SM contribution arises 
only as a second order effect. Thus, those NP contributions which do not vanish 
in the flavor SU (3) limit must determine the lifetime difference there, even if their 
contibutions are tiny in the individual decay amplitudes A simple calculation 
reveals that NP contribution to y can be as large as several percent in R-parity- 
violating SUSY models or as small as ^ 10"^" in the models with interactions 
mediated by charged Higgs particles^^. 

Consider, for example, models that do not vanish in the SU{3) limit. The two 
most common scenarios involve models whose low energy | AC| = 1 effective hamilto- 
nian involves (V-A)g)(y-A) and (S-P)®(S-|-P) couplings. For instance, SUSY mod- 
els without R-parityl2l would fit the bill. In these models, there are flavor-changing 
interactions of sleptons that can be obtained from the lagrangian 

^ KjkLiQjDl , (12) 

The slepton-mediated interaction is not suppressed in the flavor SU(3) limit and 
leads to 

C'A 

Ml 



{C2~2Ci){Q') + {Ci-2C2){Q') , (13) 



where C = -GFml/{<oV2-KMDTD), Mj is a slepton mass, A is given by A = Xsd — 
A {Xdd - Xss) - y? (\ds + A^d) with Xsd = A-12A-21 < 1 x IQ-^, A^^ = y^;,Y\Ki\ < 
5 X 10-5, \dd = < 5 X 10-^ \ds = AJnA-22 < 5 x 10^^ (see [B]), and (Q') 

is 

(QO - (D°|u,7^PiC, ?Ij7'^PflC,|75") . (14) 

Operators with a tilde are obtained by swapping color indices in the charm quark 
operators. Using factorization to estimate matrix elements of the above operators 
and taking for definiteness Mj— 100 GeV, we arrive at j/^j ~ — 3.7%P=21. 

On the other hand, there are also several reasons that some NP models vanish 
in the flavor SU{3) limit. First, the structure of the NP interaction might simply 



2, 2008 8:0 WSPC/INSTRUCTION FILE petrov 



CHARM MIXING IN THE STANDARD MODEL AND BEYOND 7 

mimic the one of the SM. Effects Hke that can occur in some models with extra space 
dimensions. Second, the chiral structure of a low-energy effective lagrangian in a 
particular NP model could be such that the leading, mass-independent contribution 
vanishes exactly, as in a left-right model (LRM). Third, the NP coupling might 
explicitly depend on the quark mass, as in a model with multiple Higgs doublets. 
However, most of these models feature second order S'C/(3)-breaking already at 
leading order in the l/rric expansion. This should be contrasted with the SM, where 
the leading order is suppressed by six powers of rris and the second order only 
appears as a 1/m^-order correction. 

For instance, LR M p rovide new tree-level contributions mediated by right- 
handed {W'^^^) bosons"'^^'. The relevant effective lagrangian is 

£lR = -^V^^^ Ua,^rPBdb,^W|,''^ +h.C. , (15) 

where V^^-* are the coefficients of the right-handed CKM matrix. Since current 
experimental limits allow Vl^(^) masses as low as a a sizable contribution to 

y is quite possible. We obtain 

2/lr = -ClrVWvW* [Ci(Q')+C2(Q')] , (16) 

where Clr = XCf^GFrnlxs/iTrMuTu), G?VV2 = g^/SM^^, Ci,2 are the SM 
Wilson coefficients and the operators appearing in Eq. (fTBl) are given in Eq. (|14p . 
Using [4], we obtain numerical values for two possible realizations: (i) 'Manifest LR' 
(V(i^) = V'^)) gives 2/LR = -4.8 • 10"^ with Mw,,, = 1.6 TeV and (ii) 'Nonmanifest 
LR' (V^^^^^ ~ 1) gives yLR = -8.8 • IQ-^ with Mw^ = 0.8 ^e^^^. 

As one can wee, small NP contributions to |AC| = 1 processes produce substan- 
tial effects in the D^Tf lifetime difference, even if such contributions are currently 
undetectable in the experimental analyses of charmed meson decays. Coupled with 
a known difhculty in computing SM contributions to D-meson decay amplitudes, it 
might be advantageous to use experimental constraints on y in order to test various 
NP scenarios due to better theoretical control over the NP contribution and SU(3) 
suppression of the SM amplitude. 

Nevif Phyiscs and CP-violation in D^D mixing. Another possible man- 
ifestation of new physics interactions in the charm system is associated with the 
observation of (large) CP-violation. This is due to the fact that all quarks that 
build up the hadronic states in weak decays of charm mesons belong to the first two 
generations. Since 2x2 Cabbibo quark mixing matrix is real, no CP-violation is 
possible in the dominant tree-level diagrams which describe the decay amplitudes. 
CP-violating amplitudes can be introduced in the Standard Model by including 
penguin or box operators induced by virtual 6-quarks. However, their contributions 
are strongly suppressed by the small combination of CKM matrix elements VcbV*^^. 
It is thus widely believed that the observation of (large) CP violat ion in charm 
decays or mixing would be an unambiguous sign for new 

physicPini. 
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5. Summary 

Charm physics provides new and unique opportunities for indirect searches for New 
Physics. NP can affect charm mixing in a variety of ways, mainly affecting both 
X and y, as well as providing CP-violating asymmetries. Expected large statistical 
samples of charm data will allow new sensitive measurements of charm mixing and 
CP-violating parameters. 
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